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Abstract: We explicitly construct all supersymmetric flux vacua of a particular 
Calabi-Yau compactification of type IIB string theory for a small number of flux 
carrying cycles and a given D3-brane tadpole. The analysis is performed in the large 
complex structure region by using the polynomial homotopy continuation method, 
which allows to find all stationary points of the polynomial equations that charac- 
terize the supersymmetric vacuum solutions. The number of vacua as a function of 
the D3 tadpole is in agreement with statistical studies in the literature. We calcu- 
late the available tuning of the cosmological constant from fluxes and extrapolate to 
scenarios with a larger number of flux carrying cycles. We also verify the range of 
scales for the moduli and gravitino masses recently found for a single explicit flux 
choice giving a Kahler uplifted de Sitter vacuum in the same construction. 
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1 Introduction 

The question of the vacuum states of string theory has been an outstanding prob- 
lem in the field since its founding days. Those vacua which could describe our 
world are of particular interest. Apart from a standard model sector including the 
corresponding gauge groups and matter in the right representations of these gauge 
groups, such vacua should have a positive vacuum energy that is extremely small 
in terms of the natural scale of gravity, the Planck / string scale [1-3]. In addi- 
tion, one may wish for 4D = 1 supersymmetry for its power to keep control over 
quantum corrections to the theory above the supersymmetry breaking / electroweak 
scale. There has been remarkable progress in the construction of such de Sitter (dS) 
vacua with stabilized geometric moduli in the past years [4-6]. The moduli are 4D 
massless scalar fields parametrizing the geometric deformation modes of the compact 
six-dimensional spaces all viable Kaluza-Klein type vacua of string theory need to 
describe our effectively four-dimensional reality. 

A corner of the string theory landscape where moduli stabilization can be ad- 
dressed very explicitly is type IIB compactified on orientifolded Calabi-Yau three- 
folds. The four-dimensional effective action of the geometric moduli is given by an 
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A/" = 1 theory of a set of chiral multiplets consisting of the axio-dilaton, h^'^ Kahler 
moduh, and h'^'^ complex structure moduh [7]. Recent years saw a progress for the 
stabihzation of the dilaton and the complex structure moduli from the use of quan- 
tized fluxes of three form field strength of the Ramond-Ramond and Neveu-Schwarz 
sector of type JIB string theory [5] . The flux stabilization procedure operates super- 
symmetrically at a high scale. The Kahler moduli are fiat directions at tree level, i.e. 
they do not obtain a scalar potential due to the no-scale structure of type IIB com- 
pactified on Calabi-Yau [8, 9]. We can use the breaking of this no-scale structure by 
non-perturbative [6] and perturbative [10, 11] effects to stabilize the Kahler moduli 
at a parametrically lower scale than the complex structure moduli. This stabilization 
produces in an AdS vacuum with unbroken supersymmetry [6], an AdS vacuum with 
spontaneously broken supersymmetry with an exponentially large volume [12, 13] or 
directly in a dS vacuum with spontaneously broken supersymmetry [14, 15]. 

If the stabilization of the geometric moduli does not directly lead to a dS vacuum, 
an additional uplifting sector has to be added. The uplifting may arise for instance by 
anti D3 branes [6], D-terms [16-18], F-terms from matter fields [19], metastable vacua 
in gauge theories [20] or dilaton dependent non-perturbative effects [21]. Concerning 
the smallness of the cosmological constant, statistical arguments show that due to the 
enormous amount of possible flux configurations there is an exponential abundance 
of isolated potential dS vacua [22-24]. For a sufficiently large number of complex 
structure moduli h"^'^ (typically (9(100)) the number of flux vacua and, in turn, also 
the number of dS vacua scaling like e'-^^^^^^'^ is large enough to produce vacuum 
energies with average spacing < 10"^^°. This enables the flux vacuum landscape 
in string theory in principle to accommodate the observed vacuum energy of our 
Universe. 

In this paper, we study the flux vacua of a particular Calabi-Yau: The degree 18 
hypersurface in a 4 complex dimensional projective space Xs = CPtn69[18] in the 
large complex structure limit. This manifold is the standard working example of both 
the large volume scenario (LVS) [12, 13] and the Kahler uplifting scenario [14, 15] 
and its geometric properties have been worked out in great detail in [25]. It has 
h^'^ = 2 Kahler moduli and /i^'^ = 272 complex structure moduli. We switch on 
flux along six three-cycles that correspond to two complex structure moduli that are 
invariant under a certain discrete symmetry that can be used to construct the mirror 
manifold [26]. For this purpose we review a known argument that a supersymmetric 
vacuum in these two complex structure moduli corresponds to a supersymmetric 
vacuum of all 272 complex structure moduli [27, 28]. 

For an explicit construction of the flux vacua we use the fact that the prepotential 
Q of the two complex structure moduli space has been worked out in [25] in the large 
complex structure limit. We apply two computational methods to find flux vacua on 
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this manifold: 

• The polynomial homotopy continuation method [29] allows us to find all station- 
ary points of the polynomial equations that characterize the supersymmetric 
vacuum solutions. The fluxes ji ^1- appear as parameters in these equations 
and are restricted by the D3 tadpole L which depends on the chosen brane 
and gauge flux configuration imposed on the manifold. Since the restriction 
is of the form X] — ^'^is method allows us to explicitly construct for the 
first time all flux vacua in the large complex structure limit that are consistent 
with a given D3 tadpole L by applying the polynomial homotopy continuation 
method at each point in flux parameter space. This method has the attractive 
feature to be highly parallelizable. 

• The minimal flux method [28] finds flux parameters that are consistent with 
a given D3 tadpole L for a given set of vacuum expectation values (VEVs) of 
the complex structure moduli. Hence it is in a sense complementary to the 
polynomial homotopy continuation method where the role of parameters and 
solutions is exchanged with respect to this method. However, it is not possible 
to find all flux vacua for a given tadpole L with this method. 

The obtained solution space of flux vacua is analyzed for several physically inter- 
esting properties. For the polynomial homotopy continuation method, we find that 
for the ~ 50, 000 flux choices contained in our maximum D3-brane tadpole L = 34 of 
our scan there are ~ 20, 000 solutions in the large complex structure limit. We find a 
preference of strongly coupled vacua gs ^ ^ and preference for values of C(10^ — 10^) 
for the flux superpotential Wq. The number of vacua is 

N^ac ^ (0.52 ± 0.04) l2-92±o.o3 ^ 

compared to ~ 0.03L^ expected form statistical analysis [24, 30]. The gravitino 
mass is typically ml^^ = C(10~'^) ■ (^)^M| and the masses of the complex structure 
moduli and the dilaton scale like C(10"^ — 10^) • {^^Mp, where V is the Volume 
of X3 in string units. 

The average spacing of the flux superpotential in our solution set can be used to 
estimate the available fine-tuning A A/A of the cosmological constant A as 

^ ^ (5.1 ± 0.3) L-(°-93±o.oo6) (hf^Vi) ^ (12) 

2 1 

where /i^g is the number of complex structure moduli with non-zero flux on the cor- 
responding three-cycles. Eq. (1.2) is obtained fit for L < 34 and h;i^ = 2. It 
can be used to estimate the available fine-tuning of the cosmological constant to for 
instance AA/A ~ 10^^^^ for L = 500 and hl^ = 40. These are typical values for 
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Calabi-Yau manifolds that are hypersurfaces in toric varieties with D7 branes and 
07 planes introduced to stabilize the Kahler moduli and break supersymmetry. 
The explicit brane and gauge flux construction in [15] allows us to answer the ques- 
tion how many of theses supersymmetric flux vacua allow an uplift to dS via Kahler 
uplifting. Depending on the available values for the one-loop determinant from gaug- 
ino condensation used to stabilize the Kahler moduli in this setup, we find that for a 
fraction of about 10~^ of all flux vacua up to a given D3-brane tadpole this mechanism 
can be applied to obtain a dS vacuum. 

For the minimal flux method, we find ~ 1000 flux vacua with L < 500 out of 
~ 10^ parameter points of our scan. This method allows us to control the region in 
Wq and moduli space where we are intending to find flux vacua. Hence, we more 
easily access the regions of weak string coupling and the large complex structure 
limit compared to the polynomial homotopy continuation method. For this much 
smaller set of flux vacua constructed with the minimal flux method, the fraction of 
Kahler uplifted dS minima is about 10%. This is a considerably higher fraction of 
vacua compared to the polynomial homotopy continuation method which is due to 
the fact that the minimal flux method naturally finds values for Wq in a region where 
Kahler uplifting is applicable. 

Our results are complementary to statistical analysis by [24, 30].^ The uniform 
distribution of physical quantities as for instance the gravitino mass and the vacuum 
energy density in the landscape has recently been questioned in general [34-36] and 
in the context of Kahler uplifting [37-39]. Hence, our results present an important 
check of the general results found in [24, 30] on a very realistic examples X3, especially 
since we are able to construct the complete solution space of flux vacua for a given 
tadpole L. 

In section 2, we review the effective 4d, M = 1 description of the complex struc- 
ture moduli and the dilaton as well as the reduction of the full moduli space to 
two complex structure moduli. The scans for flux vacua with the polynomial homo- 
topy continuation method and the minimal flux method are presented in section 3 
respectively section 4. We conclude in section 5. 

2 The complex structure of CP^^iegt^^] 

The effective 4d, M = 1 description of the moduli space of the K^'^ = 272 complex 
structure moduli of is given by the Kahler potential K and superpotential Wq of 

^For explicitly constructed vacua on two different two parameter models in the vicinity of the 
Landau-Ginzburg respectively conifold point see [31]. For a study of flux vacua of X3 in the context 
of accidental inflation [32] see [33]. 
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the theory.^ We choose a symplectic basis {A^, B^} for the 63 = 2h'^'^ + 2 three-cycles 



a'= a''APa = Sl / /3a= / A a'' = -5^ (2.1) 

Aa JX3 JBb JXs 

where {at,/?"} are the Poincare dual cohomology elements to the three-cycles and 
a,b = 0,..,h^'\ 

Having chosen a symplectic basis for the three-cycles, this defines a choice of 
coordinates Ua on complex structure moduh space via the period integrals over the 
holomorphic three-form Q via 



u;a= g, = (2.2) 

Note, that there are h?'^ + 1 coordinates wq, ••, ^^2,1 even though there are only 
K^'^ complex structure moduli. This is because uq refers to the normalization of 
the holomorphic three-form Vt. The complex structure moduli can be defined via 
Ua = I'a + = i^a/'-^o for i = a,..,h'^'^. The period vector Il{cOa) = {Qb,<^a) is 
inherited from a holomorphic function Q{u}a) of degree two in the Ua known as the 
prepotential via Qb = dbQ of the underlying Af = 2 Calabi-Yau compactification. 

The Kahler potential of the complex structure moduli Ua and the dilaton r = a+i s 
can then be written as 

i^cs = - log (^-^ niUa) A - log (-z(r - f )) , 

= - log {l{UaQa - i^aQa)) - log {-i{t - f )) , (2-3) 

= - log (-znt ■ S • n) - log (-z(r - f )) , 
where in the second line of eq. (2.3) we have introduced the symplectic matrix 

1 



-1 



(2.4) 



and used the intersection formula 

/ XAY = J2([ x[ Y- [ y[ x), (2.5) 

JXs a=0 "^-S" JAa JB'^ J 

for general three-forms X and Y . 

The Gukov-Vafa-Witten flux superpotential is determined by the RR and NS flux 
F3 and ifa via [43] 

Wo = ^ [ {Fs-rH3)AQ{Ua). (2.6) 



^For recent reviews see [40-42]. 



- 5 - 



Due to the quantization of the three-form flux 



^ ^ F3 = /i,GZ, jAn I Fs = f2aez. 



(2.7) 



(2vr)V7^^ ^ '-''^ ' (27r)V 

— ^ / H3 = hiaez, \ ^ I i/g = /i2a e Z , 

(27r)^a' (27r)^«' 

eq. (2.6) can be written as 

Wo = 2n [(/i, - T h,)g, - U2a - r haWa] , (2.8) 

where we have set a' = 1 and used again eq. (2.5) and the definition of the periods 
eq. (2.2). 

The D3-tadpole induced by turning on RR and NS flux is given by 

L = , ,i I H3AF3 = h-E- f = hj2- h2fi . (2.9) 
(27r)4(a')2 J^^ ■> '•'^ ^ ^ 

The M = 1 supergravity scalar potential is given as 

V = (^K'^^DaWD^ -3\W\'^^ , (2.10) 
where K = K^s + Ky^ and 

irk = -21ogV, (2.11) 

is the Kahler potential of the Kahler moduli up to corrections in a' and Qs with V 
the volume of the Calabi-Yau X^. The indices a and /3 in eq. (2.10) run over the 
dilaton, the /i^'^ complex structure moduli and the h^'^ Kahler moduli. At tree-level, 
eq. (2.10) obeys a no-scale structure [8, 9] in the Kahler sector: 

K'WiW'D]W = 3\W\^ , for 2,j = l,..,/i^'S (2.12) 

such that 



^ = K'^D^WDdW , (2. 13) 

where the indices c and d run over the moduli r and Ua- The no-scale structure 
eq. (2.12) is broken by a' corrections [10] and string loop corrections [11] in as 
well as non-perturbative corrections to the superpotential and K. However, these 
corrections are parametrically small in every moduli stabilization scenario where the 
overall volume of X3 is large. Hence, the scalar potential for the dilaton and complex 
structure moduli eq. (2.13) is positive semi-definite in this limit and a supersymmetric 
extremum given by a solution to the system of equations 

D^W = Q and Du^W = Q , for a = 1, .., /i^'S (2.14) 
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will always be a minimum, i.e. all eigenvalues of the second derivative matrix Vat are 
positive [13]. 

Note that due to the appearance of the symplectic matrix, the tadpole eq. (2.9) 
is at first not positive definite. However, as has been discussed in [5, 24], imposing 
the supersymmetry conditions DaW = results in G3 = -F3 — r being imaginary 
self-dual (ISD). Since the ISD component of G3 always results in positive semi- 
definite contributions to the tadpole while the anti-ISD component of G3 always 
yields negative semi-definite contributions, a supersymmetric point always has L > 0. 
This can be seen nicely if the ISD condition is displayed as [44] 

HsHs = -{Fs-aHs), (2.15) 

i.e. only + 1 of the original 2h'^'^+2 fluxes are independent once the ISD condition 
is invoked and 

L~ /" H^AF^r^ f H^A^H^r^ j ^g\H^\^ > Q , (2.16) 

J X3 J X2, J 

where we have used eq. (2.15). 

The type IIB ten dimensional effective supergravity Lagrangian is invariant under 
SL{2, Z) transformations 



ar + b 
CT + d 



T — )• ^ — ; with a,b,c,dG Z and ad — be = 1 , 



H3\ d c\ fH, 
fJ [b a '[Fs 



(2.17) 



which implies 



Gs^-^. (2.18) 

CT + d 

As is easily verified, these transformations also leave the D3 tadpole eq. (2.9) invari- 
ant. When determining the solution space of flux vacua of X3 we have to make sure 
to consider only inequivalent vacua under the transformations eq. (2.17). 



2.1 Effective reduction of the moduli space 

Consider the two parameter 0-family of threefolds CP^;^]^g9[18] given by the van- 
ishing of the polynomials 

+ + ^3*^ + X4 + X5 — l^ipXiX2X'iX4^x^ — ?>(l)x\x\x^ , (2.19) 

i.e. all except two of the 272 complex structure moduli which correspond to monomi- 
als in the general degree 18 Calabi-Yau hypersurface equation have been set to zero. 
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As was discussed in [25], eq. (2.19) is invariant under a global F = Zg x Zig sym- 
metry. This symmetry is used in the Greene-Plesser construction [26] to construct 
the mirror Calabi-Yau which in this case has h^'^ = 272 and h'^'^ = 2. Furthermore, 
the moduli ip and in eq. (2.19) describe the two complex structure moduli of this 
mirror manifold. As was pointed out in [45], the periods of the mirror agree with 
those of CP^;^]^g9[18] at the F symmetric point. Also, [45] shows that the complete set 
of h"^'^ complex structure moduli can be divided into a F-invariant subspace and its 
complement. The moduli with trivial transformation are exactly those that do not 
vanish at the F symmetric point, in this case ip and (p. 

To make use of the agreement of the prepotential for the complex structure sector 
of CP^]^ig9[18] and its mirror in the large complex structure limit, it is useful to 
introduce the complex coordinates Ui and U2 that are related to ip and (p as [25] 

Xi = -— (1 + 312gi + 2q2 + 10260^,^ - 540gig2 - ql 

- 901120^1^ + 120420g^g2 + 20g^ + . . . ) , 
^ (2.20) 

X2 = (1 + 180gi - 6q2 + llGlOg^ + 180gig2 + 27^2^ 

+ 514680^1^ - 150120g^g2 - 5040gig^ - 164g^ + . . . ) , 
up to third order in the qj = e^'^*^-' with the large complex structure coordinates 

^^^(18^ and X2 = (30)^ (2.21) 

6(p 

The large complex structure limit corresponds to Xj — )■ 00 which is equivalent to 
lm{Uj) — )■ 00 as can be seen from eq. (2.20). 

There are two conifold singularities given by the equations [25] 
CFi : (26244^^*^ + <P? = I ^ / ^ 



CFo: 6^ = 1 ^ 



X2 


(- 


27 


V432 


X2 






= 1. 


27 





(2.22) 



Let us come back to the problem of finding supersymmetric extrema by solving 
eq. (2.14). As was noted in [28, 46], to find an extremum it is sufficient to turn on 
fluxes only along the six F-invariant three-cycles, i.e. 

/ = (/ii,/i2>/i3'/2i>/22,/23,0, ...,0) and h = {hi^,hi^,hi^,h2,,h22,h2,,0, ...,0) , 

(2.23) 

having set to zero all the components along the &3 — 6 non-invariant three-cycles. 
It is then possible to achieve DaW = for all 272 complex structure moduli,'^ and 



^Note that orientifolding will project out some of the 272 complex structure moduli. Since the 
exact number of projected out directions depends on the position of the 0-plane we stick to the 
upper bound of 272 for a general treatment. 
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hence to find a minimum of the positive definite tree-level no-scale scalar potential 
eq. (2.13). This is due to the fact that, for this F invariant flux, the symmetry F is 
realized at the level of the four- dimensional effective action. Note that the restriction 
to flux on the F invariant cycles is purely for simplicity, as the analysis of the complete 
272 dimensional complex structure moduli space is practically extremely challenging. 

Let us explain more detailed why the flux vector in eq. (2.23) generically provides 
a stable minimum of all 272 complex structure moduli [28, 46]. We flrst consider 
D^^Wq = 0, where for a = 3, . . . , 272 denote the non-trivially transforming mod- 
uli under F = Zg x Zig. In the large complex structure limit, the prepotential Q 
is a polynomial function of all h'^'^ complex structure moduli that has to transform 
trivially under F, since if it would not, F could be used to flx the non-trivially trans- 
forming moduli.^ Hence, the non-trivially transforming Ua have to appear at least 
quadratic in Q in order to represent a F invariant contribution to Q. This informa- 
tion, together with having switched on flux only along the F invariant directions, is 
sufficient to show 

W^ fj^ = Kfj^ = at Ua = for a = 3, . . . , 272 , (2.24) 

since Wq is a polynomial function which is at least linear in the see eq. (2.8) 
and Kjj is a rational function which is at least linear in the numerator in the Ua, 
see eq. (2.3). Hence, D^jWo = W^ fj^ + K^Wq = at Ua = for a = 3, ... , 272. 
This reduces the full set of conditions DaW = Va to the three equations 

DiW\^^^Q = for J = r, Ui, U2 . (2.25) 

This is equivalent to set Ua = from the beginning and study the stabilization 
problem for the reduced case with two complex structure moduli, as we do in the 
following. 

In [25], the prepotential Q for the two complex structure moduli Ui = Ui/uq and 
U2 = was derived via mirror symmetry in the large complex structure limit to 

be 

^(wo,wi,w2) = ^u;o+ — ^ — ^^2r^^ 6^;^ ' ^ ' 

with = 2{2T!i)'i — — 1-30843 z determined by the Euler number x of the Calabi-Yau. 
Eq. (2.26) receives instanton corrections which are given as 

1 / 1215o^ 
QinstXqi, 52) = 77— 540gi + — ^ + 560g? + 3^2 - 1080gig2 + 143370g2g2 

_15£l + 2rOO,,,,^ + ?!M + ...V (2.27) 
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'^Q completely determines the moduli space of the (before orientifolding) M = 2 moduli space. 
If it would not be invariant the complex structure moduli space would have been reduced, i.e. some 
flat directions lifted but this does not happen just because there exists a T symmetric point. 
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with Qa = exp {271 iUa) and we have set Uq = 1. The dots in eq. (2.27) denote 
higher powers in the which are suppressed in the large complex structure limit 
Ua = lm{Ua) > 1. We define the large complex structure limit via 



/inst. 



< ^LCS : 



540e- 



-2'rrui 



< ^LCS 



and 



3e 



-2lTU2 



(2.28) 



\g\ {27ry\g\ - (2vr)3|e?| 

for small e^cs- The two last condition in eq. (2.28) are imposed to ensure that there 
are no cancellations between the terms in ^inst., i-e. the leading correction in e"^'^"" 
is actually small. Furthermore, to have a valid description of the complex structure 
moduli by g we have to ensure that we are not in the vicinity of the conifold points 
eq. (2.22), i.e. 

3 



X2 

27 



432 



+ 1 



1 



and 



X2 

27 



1 



(2.29) 



with small ecF- 



3 The polynomial homotopy continuation method 

We want to solve the non-linear eqs. (2.25) derived from the prepotential eq. (2.26) 
for the 6 real variables Xj = Mi, M2, s, ui, 1^2 and a. The parameters of these equations 
are the 12 fiuxes /i, /2, hi and /12 in eq. (2.23). Though systems of non-linear equa- 
tions are extremely difficult to solve in general, if the non-linearity in the system is 
polynomial-like, then the recently developed algebraic geometry methods can rescue 
the situation. In particular, we use the so-called numerical polynomial homotopy 
continuation (NPHC) method [29] which finds all the solutions of the given system 
of polynomial equations. This method has been used in various problems in particle 
theory and statistical mechanics in Refs. [47-57]. 

3.1 The algorithm 

Here we briefiy explain the NPHC method: for a system of polynomial equations, 
P{x) = 0, where P{x) = {pi{x), . . . ,Pm{x))'^ and ) , which is known 

to have isolated solutions, the Classical Bezout Theorem asserts that for generic 
values of coefficients, the maximum number of solutions in C" is Hiii'^j- Here, di 
is the degree of the ith. polynomial. This bound, the classical Bezout hound (CBB), 
is exact for generic values [29, 58] for details. 

Based on the CBB, a homotopy can be constructed as 

H{x,t) = -i{l-t)Q{x)+t P{x), (3.1) 

where 7 is a generic complex number and t G [0, 1). Q{x) = (gi(x), . . . , qm{x))'^ is a 
system of polynomial equations with the following properties: 
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1. the solutions of Q{x) = H{x, 0) = are known or can be easily obtained. Q{x) 
is called the start system and the solutions are called the start solutions, 



2. the number of solutions of Q{x) = H{x, 0) = is equal to the CBB for P{x) = 



3. the solution set of H{x,t) = for < t < 1 consists of a finite number of 
smooth paths, called homotopy paths, each parameterized by t G [0, 1), and 

4. every isolated solution of H{x,l) = P{x) = can be reached by some path 
originating at a solution of H{x, 0) = Q{x) = 0. 

The start system Q{x) = can for example be taken to be 



where di is the degree of the i^^ polynomial of the original system P{x) = 0. Eq. (3.2) 
is easy to solve and guarantees that the total number of start solutions is YYiLi 
all of which are non-singular. 

We can then track all the paths corresponding to each solution of Q{x) = from 
t = to t = 1. The paths which reach P{x) = = H{x, 1) are the solutions of 
P{x) = 0. By implementing an efficient path tracker algorithm, all isolated solutions 
of a system of multivariate polynomials system can be obtained because it is shown 
[29] that for a generic 7, there are no singularities (i.e., paths do not cross each other) 
for t G [0, 1). In this respect, the NPHC method has a great advantage over all other 
known methods for finding stationary points. 

There are several sophisticated numerical packages well-equipped with path track- 
ers such as Bertini[59], PHCpack [60], PHoM [61] and HOM4PS2 [58, 62]. We mainly 
use Bertini to get the results in this paper. 

We mean by a solution a set of values of variables which satisfies the eqs. (2.25) 
with tolerance 10~^°. All the solutions come with real and imaginary parts. A 
solution is a real solution if the imaginary part of each of the variables is less than 
or equal to the tolerance 10~^ (below which the number of real solutions does not 
change, i.e., it is robust for the problem at hand). All these solutions can be further 
refined to an arbitrary precision. 

The advantages of the homotopy based on the CBB are (1) the CBB is easy to 
compute, and (2) the start system based on the CBB can be solved quickly. The 
drawback of it is that the CBB does not take the sparsity of the system into account: 
systems arising in practice have far fewer solutions than the CBB, so a large portion 
of the computational effort is wasted. 




1\ 



Q{x) 



(3.2) 
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Hence, one can also use homotopies based on tighter upper bounds. For exam- 
ple, one can compute the so-called 2-Homogeneous Bezout Bound or the Bernstein- 
Khovanskii-Kushnirenko bound [63-65] which are tighter upper bounds. These two 
bounds were explained in Ref. [47, 54]. We note that, as with the CBB, the 2HomBB 
and BKK bound are also generically sharp with respect to the family of polynomial 
systems under consideration. 

There is yet another, rather more practical, way of solving a parametric system 
which is called Cheater's homotopy [66, 67]: let us say we want to solve a parametric 
system, f{q; x) = where x are variables and q are parameters, in our case the fluxes. 
Now, it can be shown [66, 67] that the maximum number of complex solutions at 
any parameter point is the number of solutions at a generic parametric point. So our 
strategy is first to solve the system at a generic parameter point and then using the 
solutions at this point as the start solutions for the systems at all other parameter- 
points. This homotopy is called cheater's homotopy. A recently developed software, 
based on Bertini, called Paramotopy [68], precisely uses cheater's homotopy and 
goes over a huge number of parameter points in parallel. The package is publicly 
not available yet though the respective research group has kindly given access to the 
code for the purpose of the computation in this paper. We will publish the details 
on the cheater's homotopy and the package Paramotopy elsewhere [69]. 

3.2 The scan 

We define a set of flux parameters on which we apply the algorithm described in the 
previous section 3.1. Since we are only interested in supersymmetric flux vacua, we 
can make use of the ISD condition eq. (2.15) and define a flux configuration via 

with hi,h2 G Z^. Note that since we have two complex structure moduli we have 
initially 2-2 + 2 = 6 flux parameters for both and -F3 but the ISD condition 
eq. (2.15) reduces this to the six parameters given in eq. (3.3). Furthermore, the D3 
tadpole eq. (2.9) becomes manifestly positive semi-definite, i.e. 

L = hl + hl. (3.4) 

To scan efficiently, we apply the paramotopy algorithm only to SL{2, Z) inequiv- 
alent flux configurations. Note that a configuration of the form eq. (3.3) transforms 
as 



- 12 - 



under SL{2,Z) transformations, eq. (2.17). For general a,b,c and d, F3 in eq. (3.5) 
will not be of the form in eq. (3.3) but only for the 4 cases 

a = ±1, 6 = 0, c = 0, d = ±1 and a = 0, b = ±1, c = T^, d = . (3.6) 

This corresponds to the SL{2, Z) equivalent flux configurations 

Note that the two transformations on the LHS of eq. (3.6) do not transform the 
dilaton r' = r while the two transformations on the RHS act as r' = — 1/t. 

The number of SL{2, Z) inequivalent flux configurations in a spherical region 
defined by a spherical constraint eq. (3.4) can be estimated as 7r^/(4 r[4])(v^)^, 
using the formula for the volume of the n-sphere Vn{r) = '7r"/^/r(n/2 + l)r". The 
factor 1/4 accounts for the 4 equivalent configurations in eq. (3.7). If we had switched 
on more flux n > 6 the number of lattice points grows very rapidly ~ L"/^. 

For our scan, we choose L = 34 such that we scan over 52,329 parameter points 
(the above estimate yields 55,391). On the FermiLab cluster using 100 nodes each 
with 32 cores (each core with 2.0 GHz cloak speed), the calculation time in total is 
around 75, 000 hours, with 60 — 100 minutes per parameter point. 



3.3 Distribution of parameters 

In this section, we want to discuss the distribution of the following parameters as 
results of the scan defined in the previous section 3.2: 

• Ui and U2, to see how many points reach the large complex structure limit as 
defined in eq. (2.28). 

• r, to identify regions of weak respectively strong coupling. 

• The number of solutions for a given D3 tadpole L. 

• Wq, the flux superpotential. 

• The masses of the moduli and the gravitino mass "^3/2- 

• The available fine-tuning AA of the cosmological constant A. 

• The amount of flux vacua for which a dS vacuum can be constructed via Kahler 
uplifting in 3.4. 
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-10 5 5 10 

Ui 

Figure 1: Distribution of ui and U2 for the complete set of solutions to eqs. (2.25) 
(gray) and for the physical solutions fulfilling the criterion for the validity of the large 
complex structure limit (blue for eLcs = 10~^ and black for eLcs = 10"^), eq. (2.28). 



For the 52,329 parameter points, we find a total of 500,865 solutions to the 
eqs. (2.25). This corresponds to an average of 9.5 solutions per parameter point. 
For 1,270 parameter points we do not find a solution. Many of the solutions are 
unphysical and hence have to be sorted out: A subset of 271,825 fulfill the criterion 
of a physical string coupling Qs > and only a subset of 24,882 respective 15,392 is 
in accordance with the large complex structure criterion eq. (2.28) for eics = 10~^ 
respective eLcs = 10"^. Of these none have to be sorted out because they are in 
the vicinity of the conifold singularities eq. (2.29) having chosen ecF = 10~^. Due 
to the strong suppression of the large complex structure limit in the general solution 
space of eq. (2.25) the minimal flux method has the advantage of directly searching 
for solutions in this region. 
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For the distribution of the dilaton, we can 
use SL{2, Z), to transform each solution to 
the fundamental domain 



and 



Irl > 1 



via the successive transformations 



T + b, G' 



0, d 



i.e. a = 1, b E Z, c 

-1/r, G', = Gs/t 



1 and 



r 



(3.8) 



(3.9) 



(3.10) 




i.e. a = 0, b = —1, c = 1, d = 0. 

We show the distribution of the obtained 

values for t = a + is in Figure 2. We see 

that the the strongly coupled region s = 

I/qs ~ 1 is preferred and large values of 

s > 10 are obtained for a fraction of 5%. t,. n t^- . -i r r 1.1 

l^igure 2: Distribution 01 r tor the 

paramotopy scan with eLcs = 10~^. 

The number of vacua of X3 in the large complex structure limit for a given D3 
tadpole L was estimated in [28] as ^ 

(27rL)3 



det(-7^ - 1 



UJ 



(3.11) 



with Kahler form u and the curvature two-form TZ of the moduli space. The integral 
in eq. (3.11) was estimated in [28] to be be 1/1296, using the F symmetry of the 
moduli space such that 

N^ac^Omi^. (3.12) 

Since paramotopy allows us to find all solutions for a given flux configuration we 
can not only check the L dependence of eq. (3.11) but also the normalization. This 
depends on the value chosen for eLcs: i-e- a greater eics will yield a larger normal- 
ization factor. Fitting the number of solutions with < L in the large complex 
structure limit to the tadpole L we find 



"vac — 

AT ~ 

"vac — 



(0.52 ± 0.04) L2-92±o.o3 
(0.88 ± 0.06) L2-9i±o o3 



for 
for 



^LCS 
(-LCS 



10" 

10" 



(3.13) 
(3.14) 



^Note that Nyac ^ in [28] which is due to the fact that 12 independent fluxes have been 
switched on while we effectively switch on 6 independent fluxes, see eq. (3.3). 
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Figure 3: The number of vacua N^ac with < L (left) and the logarithmic dis- 
tribution of the flux superpotential Wo (right) in the large complex structure limit 
with eics = 10"^. 



The dependence of N^ac on L and the fit in eq. (3.13) are shown in Figure 3. Consid- 
ering the very general arguments that are used to derive the estimate eq. (3.11), the 
agreement within an order of magnitude with the factual number of vacua strongly 
confirms the statistical analysis of [24, 30]. In the following, we set eics = 10~^. 

The distribution of the flux superpotential is shown in Figure 3. We find that 
for most vacua O{10^ — 10^) values are preferred. To calculate the masses of the 
moduli we have to know the value of the volume V of X3 which enters via the Kahler 
potential of the Kahler moduli given in eq. (2.11). Note that we have not specified 
the stabilization mechanism for the Kahler moduli and hence have no information 
about the value of V. For the KKLT and Kahler uplifting scenarios the volume is 
typically stabilized at O{10'^ - 10^) while for the LVS it is O{10^ - 10^^). Hence, we 
can only calculate the physical masses m up to factors of V~^, i.e. 

m = ^, (3.15) 

where rrics is the mass calculated from the effective theory of the complex structure 
moduli only, i.e.i^ = Kcs and W = Wq. 

The distribution of the physical moduli masses in terms of m^^, i.e. the 
eigenvalues of the Hessian dadhV of the no-scale potential eq. (2.13) for a,b = 
Ui,U2,s,h'i,h'2,<J is shown in Figure 4 as well as the gravitino mass rriy^ terms 
of the quantity 

ml,s/2^ml/2V' = e^-\Wo\\ (3.16) 

This quantity m^^ governs the scale of the typical AdS cosmological constant 
induced by the flux superpotential ignoring the contributions from the Kahler moduli 
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Figure 4: Masses 



Af2 



^3/2 



for a 



V=100 



of the moduli, and the gravitino mass -Jr- 
v=ioo ™p 

Calabi-Yau volume V = 100. For different Calabi-Yau volumes the masses scale as 
(^) for the moduli Mi, M2, s, i^i, z/2, a (left) and for the gravitino mass as 

v=ioo 



Ml 



''3/2 



M2 



(^)^ (right). The left plot includes all eigenvalues of the Hessian, i.e. 6 

v=ioo 

values per parameter point. 



sector. 



The distribution of rn^i^ is peaking at {niy^) = 3.5 x 10^^ ■ (^)^ with standard 
deviation 3x 10~^-(^)^. The complex structure moduli and the dilaton are stabilized 
at ~ O{10~^ — 10^) (^)^. These ranges for the moduli and gravitino masses are 
compatible with the values obtained for a single explicit flux choice in the same 
construction [14, 15]. 

The AdS respective dS cosmological constant before tuning is up to OiV) factors 
estimated to be 

A~^ = %^. (3.17) 

in the LVS respective Kahler uplifting scenarios. In particular, the tunability of 
"^cs,3/2 by three-form flux directly translates into the tunability of the cosmological 
constant via 

A "^cs,3/2 

Note that the RHS of eq. (3.18) is independent of the volume V, i.e. fine tuning of 
^^08,3/2 only has a tiny effect on the VEVs of the Kahler moduli. 
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Since the polynomial homotopy continuation method allows us to calculate all 
solutions for a given tadpole L we can estimate AA/A by determining the average 
spacing Ames, 3/2 for all values of mcs,3/2 that are to be found in a cr-interval around 
(^cs,3/2)- Since the number of vacua is given as a power-law in L, with the expo- 

2 1 

nent linear in the number of flux carrying complex structure moduli h^^ , we expect 



Ames, 3/2/ 



m, 



cs,3/2 



to be of the form 
Am, 



■cs,3/2 



c 



m, 



'cs,3/2 



(3.19) 



with C, a > 0. We can determine these parameters by fitting the LHS of eq. (3.19) 
as a function of L for = 2. Choosing a 3-a interval ^ around (mcs,3/2) we find 
the available tuning for the cosmological constant to be 

(5.1±0.3)L-(o-93±o.oo6)(/^,^^.^ 



A 



^1) 



(3.20) 



where we have included the statistical errors of the fit parameters C and a. 



Let us assume that eq. (3.20) is valid and (mcs,3/2) ~ C(10) also for larger values 
of L ~ O{10^) and larger values of /i^g ~ O{10^ — 10^). Then, we can extrapolate 
the values of the cosmological constant eq. (3.17) and its tunability to more realistic 
scenarios, see Table 1. 



^2,1 
^cff 


L 


AA/A 


2 


34 


7 • 10"^ ± 5 ■ 10"^ 


2 


500 


5 ■ 10-^ ±4 ■ 10-^ 


40 


34 


5 -10-^^ ±4 -10-^^ 


40 


500 


2 ■ 10-100 ± 2 ■ 10-100 



Table 1: The tunability AA/A of the cosmological constant for different values of 
and L with statistical errors propagated from eq. (3.20). The untuned values 
of the cosmological constant are estimated via eq. (3.17) to be (9(10-"^ — 10 
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m 



units of Mp for V of O{10^ — 10^). The first row of this table is directly calculated 
from our dataset while the last three entries are obtained as an extrapolation via 
eq. (3.20). 



To tune the cosmological constant to the accuracy given in Table 1, one has 
to make the assumption that every supersymmetric flux vacuum has no tachyonic 

^There is only a week dependence on the width of the intervaL For 5-ct the difference in 
^™cs, 3/2/"^cs, 3/2 compared to 3-a is less than 1%. 

^This assumption is reasonable when the prepotential Q is of the same structure as eq. (2.26), 
i.e. we are considering the large complex structure limit away from e.g. conifold singularities via a 
mirror construction. It may be interesting to consider such examples with h^^^ = h^^^ > 2, e.g. by 
choosing random pre- factors in a general polynomial prepotential of degree 2 in the uji . 
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directions after uplifting and stabilizing the Kahler moduli. Especially, for large 
values of h^'^ there could be strong suppressions of tachyonic free configurations [34- 
36]. In the following section, we will determine how many de Sitter vacua can be 
constructed from our dataset on X3 via the method of Kahler uplifting. 

3.4 de Sitter vacua via Kahler uplifting 

The two Kahler moduli Ti and T2 of X3 can be stabilized in a dS minimum by Kahler 
uplifting. A globally consistent D7 brane and gauge flux setup that realizes such a 
dS vacuum has been presented in [15]. The Kahler potential of Ti and T2 is given as 

3/2 



K= -21o£ 



/ _ 1 _ \ ' 1 _ 1 



(3.21) 



f) = ''Z\. (-^ - > (3-22) 



with the leading order a' correction [10] 

C(3)X 
4 V2(27r)3 

with Euler number x = 2(2 — 272) = —540. To apply the method of Kahler uplifting 
we need to balance the leading order a' correction to the Kahler potential with 
non-perturbative contributions to the superpotential. These originate from gaugino 
condensation of an S'[/(24) and 5*0(24) pure super Yang Mills from respectively 24 D7 
branes wrapping the divisors corresponding to Ti and T2. The induced superpotential 
is 

W = Wo + Aie^^^' + A2e^^^K (3.23) 

By switching on suitable gauge flux it can be shown [15] that Ai, A2 7^ 0. The induced 
D3 tadpole by this gauge flux and the geometric contributions from the D7 branes 
is L = 96 — 104. The one-loop determinants Ai and A2 depend on the complex 
structure moduli, the dilaton and also potentially D7 brane moduli. The explicit 
dependence on these moduli is unknown, however for the purpose of Kahler moduli 
stabilization the values of Ai and A2 can be assumed be constant since complex 
structure moduh are stabilized at a higher scale. Choosing Ai = A2 = 1, it was 
found numerically in [15], that the pairs of Wq and s that are suitable to realize a 
dS vacuum with small positive tree level vacuum energy lie on the curve 

To parametrize our missing knowledge of the values of Ai and A2 we introduce 
the parameter AA and the scaling relations 

Wo^Wq- AA, Ai^ Ai- AA, ^ ■ , (3.25) 



*In this case, small refers to how small we can tune {V) by choosing numerical values for Wq 
and s to a certain decimal place and is not related to the tuning of the cosmological constant. 
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Figure 5: The number of Kahler uplifted dS vacua as a function of AA (left) and 
data points (s, |W^o|) (right). Kahler uplifting can be applied in the shaded region 
{AA = 2). 

under which the position of a minimum of the potential eq. (2.10) is invariant since 
V ^ V ■ AA"^. For a given uncertainty in the one-loop determinants AA^^ < Ai < 
AA around Ai = A2 = I we can then define the criterion 



AA 



< \Wo\ < W, 



dS/ 



■AA 



and 



s > 4 



(3.26) 



for a given data point (s, \ Wo\) to allow a dS vacuum via Kahler uplifting. 

We show the number of Kahler uplifted dS vacua depending on AA in Figure 5. 
Due to the suppression of weakly coupled vacua s ^ 1 and 0{1) values of the 
superpotential, the number of vacua that can be uplifted to dS via Kahler uplifting 
is strongly suppressed. For AA = 10, only 6, i.e. a fraction of ~ lO""^ of the total 
number of flux vacua allow such an uplifting. 

The available tuning of the cosmological constant via fluxes can be estimated 
again via eq. (3.18). The Kahler moduli stabilization yields a volume of V ~ 50 [15] 
such that the untuned cosmological constant is A ~ 5 ■ 10~^ and 

AA 



A 



0.51 ±0.37, 



(3.27) 



for AA = 10. We remind the reader, that this is calculated for L = 34 which is 
the maximal value we reach in our paramotopy scan and hence less than L = 104 
which is maximally allowed by the gauge flux and D7 brane construction realized for 
a Kahler uplifted dS vacuum in our explicit CP^^]^gg[18] example. 

To summarize this section, the polynomial homotopy continuation method allows 
us to find all flux vacua for a given D3 tadpole L. The number of these vacua is well 
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estimated by the statistical analysis of [24, 30]. We find that strongly coupled vacua 
s > 1 are preferred as well as 0(10^ — 10^) values of Wq. Our results can be used 
to estimate the tunability of the cosmological constant by fiuxes and the number of 
flux vacua that can be Kahler uplifted to a dS vacuum. 

4 The minimal flux method 

In this section we describe the method to find fiux vacua that has been first used by 
Denef , Douglas and Florea [28] . In contrast to the polynomial homotopy continuation 
method described in section 3, we fix starting values for the VEVs {Ui)fix, {U2) fix 
and (r) fix and solve for the fiux values / and h. 

4.1 The algorithm 

Due to the linear dependence of Wo on / and h, see eq. (2.8), the quantities DjWq = 
Woi + KiWq for I = T,Ui,U2 are linear in these fiux vectors. Hence, if we want to 
solve the system of equations 

{Wo,DrWo,Du,Wo,Du,Wo)=0, (4.1) 

this can be written as 

M-(/,/i) = 0, (4.2) 

with M G K^^^^ for general VEVs (f/i), (f/2), (r) G C. The dimensions of M are due 
to the fact that we have 8 real equations in eq. (4.1), and there are 12 fiux integers in 
total in / and h. Note that we have also included the condition Wq = in eq. (4.1). 
In fact, we are not interested in fiux vacua where Wq is strictly zero since none of 
the well studied moduli stabilization mechanisms KKLT [6], LVS [13] and Kahler 
uplifting [12, 14, 15, 70] apply in this situation. However, eq. (4.1) will only serve as 
a starting point and we will eventually end up with vacua where Wq 7^ and 0{1). 

For M G [R^^^2 ^^j^ere is no hope to find a solution of eq. (4.2) since the fiux param- 
eters have to be integers. However, if we neglect instanton corrections induced by 
eq. (2.27) and choose rational starting values {Ui), {U2), (r) G Q + i Q in the super- 
potential, the only transcendental number in eq. (4.2) is ^ = ^^^^-'^^3 = — 1.30843. .z. 
If we approximate by a rational number ^rat, for instance ^rat = —13/10 i, we have 
accomphshed M G Q^^^^. 

Now, we can hope to find a solution of eq. (4.2) although the entries of / and h will 
be generically be quite large for generic M, since one generally expects them to be at 
least of the order of the lowest common denominator of the entries of M. This puts 
tension on the D3 tadpole constraint eq. (2.9) since generally the geometry of the 
compactification manifold and D7 brane configuration generates L ~ 10^ — 10^. We 
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use the same algorithms [71] as the authors of [28], to generate as small as possible 
values for the entries of / and h in order to generate a not too large D3 tadpole 
L < Lmax where we choose Lmax = 500 to be the maximal value for the D3 tadpole 
that we consider. 

Since the system of equations eq. (4.2) is under determined, the solution space is 
given by all linear combinations of linearly independent vectors (/, h)i for i = 1, ..,4, 
where each (/, h)i is a solution to eq. (4.2), i.e. 

4 

(/, h)soi = ^ai if, h)i with G Z . (4.3) 

i=l 

For obvious practical reasons, we cannot consider all possible values of the Oj. Since 
L{ai{f, h)i) = afL{{f,h)i) and tadpoles of solution basis vectors smaller than ten 
are extremely rare, we can safely restrict ourselves to — 3 < < 3 in order to fulfill 

In the next step, we are looking for solutions to eq. (2.25), including instanton 
corrections eq. (2.27) to the prepotential and also setting ^ to its transcendental 
value. We insert the flux solution (/, h)soi into the equations 

{Dt-Wq, Du^Wq, Du^Wq) = (with instanton corrections) , (4.4) 

leaving the VEVs (Ui), {U2) and (r) unfixed. These are six real equations for six 
real variables and we can numerically search for a solution in the vicinity of (f/i) fi^c, 
{U2) fix and {t) fix- It has to be checked case by case if the complex structure limit 
is still valid for these perturbed solutions. The shift of the VEVs from their fixed 
values may also induce a shift in the superpotential, i.e. is not zero anymore. 
However, note that the value that Wq will take in the end is not in any way under 
our control and it has to be checked if one obtains useful values for the purpose of 
moduli stabilization. 

The above outlined algorithm can be iterated by sampling over a set of VEVs 
{Ui)fix, {U2)fix and (r)/i^. and approximate ^ values ^rat- 

4.2 The scan 

First of all, we need to define a set of rational starting values ((f/i), {U2), (t), Crai) fix 
over which the algorithm explained in section 4.1 can be iterated. We set the axionic 
components of the fixed VEVs to zero, i.e. 

Re(f/i) f,, = Re(f/2) fi. = Re(r) = . (4.5) 

Let X G [xmin,Xmax] represent ^rat and the imaginary parts of the moduli. We use 




p G [minz(gx™„), maxz(g ] , (4.6) 
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as the set of rational numbers that fills out the interval [xmin,Xmax]- niinz(x) and 
maxz(a;) define the smallest integer greater than or equal to x respectively the great- 
est integer less than or equal to x. Note that qmax determines how 'dense' the interval 
is filled with rational numbers. For example, one has 

|l-----2l (4 7) 

1 '4'3'2'3'4' J ^ ' 

for X e [1,2] and Qmax = 4. 







■^max 


Qmax 


# 


^rat 


0.9^ 


Lie 


20 


34 


Ui 


1 


4 


5 


31 


U2 


1/6 


4 


5 


40 


S 


1 


10 


10 


289 



Table 2: Starting values for ^rat and the imaginary parts of the moduli. The last 
column counts the number of elements yielding from the choice of Xmin, Xmax and 
Qmax according to eq. (4.6). 

The set of starting values given in Table 2 is chosen such that the starting values 
for the complex structure moduli are in the large complex structure limit eq. (2.28) 
and the string coupling = 1/s is always larger than one. The total number of 
points in the grid of starting values defined in Table 2 has 12,184,240 points. Using 
80 cores with 2.4 GHz of the DESY Theory Cluster this yields a total calculation 
time of approximately four weeks. 
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4.3 Results 



We find 1,698 solutions fulfilling the con- 
straint L < Lmax = 500 which is only 
0.01% of the total number of points of 
the scan. Hence, most of the time one 
can not find a flux vector whose entries 
are small enough to fulflll the D3 tadpole 
constraint. 

As in section 3.3, we can make use of the 
SL{2, Z) transformations eq. (3.9) and 
eq. (3.10) to transform every solution to 
the fundamental domain eq. (3.8). Iden- 
tifying equivalent solutions in this do- 
main, 1,374 elements of the original so- 
lution set are not related via SL{2, Z) 
symmetry and hence physically inequiv- 
alent. The distribution of r is shown in 
Figure 6. Compared to the paramotopy 
scan we more easily flnd weakly coupled 
vacua with s ^ 1 which is due to the fact 
that we have chosen the starting values 
for s accordingly, see Table 2. 




-0.4 -0.2 0.0 0.2 0.4 
cr 

Figure 6: Distribution of r for the 
minimal flux scan. 



The distribution of ui and U2 as well as the distribution of the superpotential are 
shown in Figure 7. We flnd that there are no flux vacua in the vicinity of the conifold 
singularities eq. (2.29) for ecF = 10^^. Also, all flux vacua fulflll the constraint of the 
validity of the large complex structure limit description, eq. (2.28) for e^cs = 10~^ 
which is again due to the chosen starting values deep in the large complex structure 
limit. Table 2. Since we solve for vanishing Wq in the flrst step of the algorithm 
eq. (4.1), we obtain a clustering around Wq = 0, see Figure 7. This is however not a 
general property of the complete solution space as we noted in section 3.3 but rather 
O{10 — 100) values are preferred. 

We flt the number of vacua as a function of the D3 tadpole L, flnding N^ac — 
0.02 L^-^^ which strongly deviates from the estimate eq. (3.11). However, our dataset 
of 1,374 is in no way representative for the total number of flux vacua with L = 500 
such that this deviation can be easily explained by insufficient statistics. 

Finally, we can calculate how many flux vacua allow a dS vacuum via Kahler 
uplifting along the lines of section 3.4. Since the minimal flux scan is setup such that 
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Figure 7: Distribution of Ui and U2 (left) and superpotential Wq for the minimal 
flux scan. 



the values of s and \Wo\ naturally lie in the region where Kahler uplifting can be 
applied we find a much milder suppression of these vacua compared to section 3.3, 
see Figure 8. For AA = 10, three of the 75 flux vacua with L = 104 and 135 of the 
1,374 flux vacua with L = 500 allow a dS vacuum via Kahler uplifting.^ Repeating 
the estimate of eq. (3.27) for the results of the minimal flux scan, the cosmological 
constant A ~ 10~^ can be tuned to an accuracy 



AA 
AA 



1.9 ±2.2 



~ 5 ■ 10"^ ± 10~^ 



for 

for 



L = 104 , 
L = 500 . 



(4. 



To conclude this section, the minimal flux method has the advantage that one can 
specify the region in parameter space where physically interesting flux vacua should 
be constructed. In our case this region is defined by the large complex structure limit, 
weak string coupling and 0{1) values of Wq. However, it is not possible to construct 
all flux vacua for a given D3 tadpole L which can be done using paramotopy and the 
method is rather inefficient in the sense that only 0.01% of the starting values yield 
a flux vacuum. 

^The maximum D3 tadpole of is L — 104 but due to the small amount of vacua we find for 
this tadpole we also show the results for L ~ 500. 
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3.0 




Figure 8: The number of Kahler uplifted dS vacua as a function of AA (top) and 
data points (s, \ Wo\) (bottom). Kahler uplifting can be applied in the shaded region 
{AA = 2) for the minimal flux scan. We show the results for L = 500 (left) and 
L = 104 (right). 

5 Conclusions 

In this paper we have studied the flux vacua of type IIB string theory compactified 
on CP^]^^gg[18] Calabi-Yau hypersurface, i.e. the standard working example of both 
the LVS and the Kahler uplifting scenario. We switch on flux along six three-cycles 
that correspond to two complex structure moduli that are invariant under a certain 
discrete symmetry that can be used to construct the mirror manifold. As explained 
in the main text, such a supersymmetric vacuum in these two complex structure 
moduli extends to a supersymmetric vacuum of all 272 complex structure moduli. 
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To explicitly construct flux vacua, we make use of the fact that the prepotential Q 
of the two complex structure moduli space has been worked out in the large complex 
structure limit. We apply two computational methods to find flux vacua on this 
manifold: The polynomial homotopy continuation method allows us to explicitly 
construct for the first time all flux vacua in the large complex structure limit that 
are consistent with a given D3 tadpole L by applying the polynomial homotopy 
continuation method at each point in flux parameter space. The minimal flux method 
finds flux parameters that are consistent with a given D3 tadpole L for a given set 
of vacuum expectation values (VEVs) of the complex structure moduli. 

We analyze the resulting solution space of flux vacua for several physically inter- 
esting properties. For the polynomial homotopy continuation method, we find that 
for the ~ 50, 000 parameter points of our scan there are ~ 20, 000 solutions in the 
large complex structure limit. We find a preference of strongly coupled vacua gs ^ ^ 
and preference for values of O{10^ — 10^) for the flux superpotential Wq. The number 
of vacua is 

N^ac ^ (0.52 ± 0.04) L2-92±0.03 ^ (5^^^) 

compared to ~ 0.03L^ expected form statistical analysis [24, 30]. The gravitino 
mass is typically rri^^^ = O{10~^) ■ {^)'^Mp and the masses of the complex structure 
moduh and the dilaton scale like C(10~^ — 10^) ■ {^)'^Mp. These ranges for the 
moduli and gravitino masses are compatible with the values obtained for a single 
explicit flux choice in the same construction [14, 15]. 

The average spacing of the flux superpotential in our solution set can be used to 
estimate the available fine-tuning AA/A of the cosmological constant A as 

^ - (5.1 ± 0.3) C'off'+i) , (5.2) 

which corresponds to for instance AA/A ~ 10^^^^ for L = 500 and /i^^ = 40. The 
explicit brane and gauge flux construction in [15] allows us to answer the question 
how many of theses supersymmetric flux vacua allow an uplift to dS via Kahler up- 
lifting. Depending on the available values for the one-loop determinant from gaugino 
condensation used to stabilize the Kahler moduli in this setup, we find that for a frac- 
tion of about 10~^ of all flux vacua up to a given D3-brane tadpole this mechanism 
can be applied to obtain a dS vacuum. 

For the minimal flux method, we find ~ 1000 flux vacua with L < 500 out of ~ 10^ 
parameter points of our scan. This method allows us to control the region in Wq and 
moduli space where we are intending to find flux vacua. Hence, we more easily access 
the regions of weak string coupling and the large complex structure limit compared 
to the polynomial homotopy continuation method. For the much smaller set of flux 
vacua constructed with the minimal flux method, the fraction of Kahler uplifted dS 
minima is about 10%. 
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